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Abstract 

Size and isolation of local populations are main parameters of interest when 
assessing the genetic consequences of habitat fragmentation. However, their rela- 
tive influence on the genetic erosion of local populations remains unclear. In this 
study, we first analysed how size and isolation of habitat patches influence the 
genetic variation of local populations of the Dupont's lark {Chersophilus duponti), 
an endangered songbird. An information-theoretic approach to model selection 
allowed us to address the importance of interactions between habitat variables, 
an aspect seldom considered in fragmentation studies, but which explained up to 
65% of the variance in genetic parameters. Genetic diversity and inbreeding were 
influenced by the size of local populations depending on their degree of isolation, 
and genetic differentiation was positively related to isolation. We then identified 
a minimum local population of 19 male territories and a maximum distance of 
30 km to the nearest population as thresholds from which genetic erosion 
becomes apparent. Our results alert on possibly misleading conclusions and sub- 
optimal management recommendations when only additive effects are taken into 
account and encourage the use of most explanatory but easy-to-measure variables 
for the evaluation of genetic risks in conservation programmes. 



Introduction 

Human activities are threatening a large number of species 
worldwide through habitat loss and fragmentation (Andren 
1994). Changes in land use often reduce the size of popula- 
tions and increase their isolation to limits where an 
increased susceptibility to stochastic factors may precipitate 
their extinction (Lande 1993; Hanski and Ovaskainen 
2000). Loss of genetic diversity, accumulation of genetic 
load and increased rates of inbreeding may reduce birth 
and increase death rates in small populations, thereby 
reducing fitness (Jaquiery et al. 2009). This effect - usually 
referred to as inbreeding depression - together with the loss 
of adaptive potential, has shown to significantly increase 
extinction probabilities, in both simulation (Saccheri et al. 
1998; Tanaka 2000; O'Grady et al. 2006) and empirical 
studies (Saccheri et al. 1998; Vilas et al. 2006). Due to the 
critical importance for species persistence and evolution. 



genetic diversity has been identified by the International 
Union for Conservation Nature (lUCN) as one of three 
levels of priority for global conservation of biodiversity 
(McNeely et al. 1990). Genetic monitoring should thus 
become a critical ingredient of conservation and manage- 
ment plans (Laikre et al. 2009). 

Two main microevolutionary processes influence genetic 
patterns in declining and fragmented populations: genetic 
drift and gene flow. These processes are determined by 
local population size and interpatch connectivity, 
parameters that are commonly taken into account in 
conservation planning under a metapopulation theory 
framework (Hanski 1999). While genetic drift causes 
random fluctuations of allelic frequencies and loss of 
genetic diversity through time as a function of effective 
population size, dispersal-mediated gene flow can buffer 
these effects in local populations. The outcome of these 
two processes is predicted under simplified equilibrium 
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models and pure-drift models. This is, however, not the 
case under the nonequilibrium situation imposed by the 
recent or on-going ft-agmentation of wildlife populations, 
where gene flow and local drift may vary over time and 
among patches. Nevertheless, larger and more connected 
populations often maintain higher genetic diversity and 
lower inbreeding and differentiation than smaller and more 
isolated populations (Frankham 1996; Cruzan 2001; Jac- 
quemyn et al. 2009). However, the effects of population 
size and isolation are usually considered independently 
(e.g. Prentice et al. 2006; Telles et al. 2007; Dillane et al. 
2008; Meeuwig et al. 2010), and only a few empirical stud- 
ies have investigated the joint effects of both parameters 
and their interaction on the genetic dynamics of wild pop- 
ulations ( Gonzalez- Varo et al. 2010; Lange et al. 2010; 
Wang et al. 2010). Contrary to theoretical expectations, 
these studies did not find this interaction to be a major fac- 
tor affecting genetic patterns. 

Critical thresholds in habitat fragmentation have been 
defined as an abrupt and nonlinear change in some param- 
eter across a small range of habitat loss (With and King 
1999). Despite of empirical evidence for critical thresholds 
in habitat fragmentation (e.g. Swift and Hannon 2010), 
they are rarely addressed in studies of genetic variability, 
inbreeding or differentiation (Cruzan 2001; Lowe et al. 
2005; Ezard and Travis 2006; Bruggeman et al. 2010; Lange 
et al. 2010). A better understanding of drift-gene flow 
interactions and fragmentation thresholds for genetic 
erosion is thus crucial to determine how human-induced 
habitat fragmentation can affect population viability and 
long-term evolutionary processes through their impact on 
genetic variation. Furthermore, the identification of critical 
thresholds should contribute to more efficient conservation 
in fragmented landscapes by helping to gauge the relative 
benefit of acting on population size or connectivity for the 
maximization of the short-term genetic viability of local 
populations. 

In this study, we evaluated the effects of population size, 
isolation and their interaction on the genetic patterns 
(diversity, inbreeding, relatedness and differentiation) of 
recently fragmented populations. Furthermore, we aimed 
to identify the thresholds of population size and isolation 
beyond which genetic erosion starts to accumulate to 
detectable levels. We focused on the endangered Dupont's 
lark (Chersophilus duponti), a markedly sedentary steppe 
specialist passerine suffering from habitat loss and frag- 
mentation in recent decades (Telia et al. 2005). Capture- 
mark-recapture methods have only detected reduced, 
short-distance movements between local populations albeit 
being separated only by few kilometres (Laiolo et al. 2008; 
Vogeli et al. 2008). The high isolation of occupied habitat 
remnants and the species' low dispersal propensity have led 
to a loss of genetic diversity and to an increase in genetic 
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structure (Mendez et al. 2011). Hence, the life history of 
our study species and its spatially structured habitat pro- 
vide a good scenario for studying the consequences of 
anthropogenic fragmentation on population genetics and 
viability. We focus in particular on those landscape and 
population variables that are easy to monitor in the field 
and are frequently used in management. By doing so, we 
aim to increase the relevance of this study for conservation 
practice and to facilitate the application of this approach to 
other species that are threatened by recent or ongoing frag- 
mentation. 

Materials and methods 

Study species and sampling 

The Dupont's lark is an endangered songbird whose habitat 
is highly restricted to steppe areas with natural vegetation 
in Spain and North Africa (Cramp 1988). The Spanish 
population has been confined to a series of fragments of 
variable size and degree of isolation, which collectively may 
hold as few as 2200-2700 breeding territories defended by 
males (Suarez 2010). The number of breeding pairs may be 
much smaller given a high male-biased adult sex-ratio in 
this species (Telia et al. 2004). Habitat loss and fragmenta- 
tion have been sufficient to extirpate many local popula- 
tions since the 1980s (Telia et al. 2005; Suarez 2010; Vogeli 
et al. 2010) and to noticeably alter cultural transmission, 
demographic and genetic patterns of the species at both 
local and regional scales (Laiolo and Telia 2005, 2006a, 
2007; Laiolo et al. 2008; Mendez et al. 2011; Vogeli et al. 
2011). This process has not occurred homogeneously 
across the species distribution, being especially intense in 
the peripheral topographical areas defined by Laiolo and 
Telia (2006a; Fig. 1). The suitable habitat for Dupont's lark 
in the core of the distribution [Iberian Mountains and Ebro 
Valley (EV)] has remained more connected, while habitat 
loss has been more pronounced in the periphery of the spe- 
cies' distribution in recent decades (Northern Plateau, 
Southern Spain and Southern Plateau; Laiolo and Telia 
2006a,b). In this sense, Dupont's lark populations are suf- 
fering a centripetal process of contraction and isolation fol- 
lowing recent fragmentation. 

During spring to autumn of 2002-2008, we attracted 
birds using playback and employed clap nets baited with 
meal worms to capture 506 adults of Dupont's larks in 36 
localities, covering the whole distribution range of the spe- 
cies in Spain (Fig. 1). Ninety- three per cent of sampled 
individuals were males, as determined by molecular analy- 
ses (Vogeli et al. 2007). We considered each steppe patch 
as the sampling unit, so areas with more patches, as the EV, 
have more sampling points (Fig. 1). A drop of blood was 
extracted for molecular analysis and stored in pure ethanol. 
AH birds were released at the site of capture, and handling 



©201 4 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd 7 (2014) 506-518 



507 



Management thresholds to avoid genetic erosion 



Mendez et al. 




Figure 1 Location of the sampled Dupont's larl< local populations in 
Spain. Different symbols indicate different population sizes (number of 
territorial males). Dark grey areas show the approximate distribution of 
the Dupont's lark in Spain based on Suarez 2010 (discontinuities among 
neighbouring small patches are not shown at this scale). The centroid 
indicates the mean geographical centre of all populations. Dashed ellip- 
ses show the topographical areas identified by Laiolo and Telia (2006a); 
EV, Ebro Valley; NP, Northern Plateau; SP, Southern Plateau; IM, Iberian 
Mountains; SS, Southern Spain. 



was specifically approved by all competing wildlife agencies 
and administrations (further details are provided in previ- 
ous studies: Vogeli et al. 2007). This research complied 
with the norms of the Spanish Animal Protection Regula- 
tion, RD1201/2005, about protection of animals used in 
scientific research, which conforms to European Union 
Regulation 2003/65. 

Genetic analyses 

DNA was extracted from blood samples following the pro- 
tocol described in Gemmell and Akiyama (1996) with two 
steps of chloroform. All samples were amplified with 14 
species-specific microsatellite loci (Mendez et al. 2011). 
None of the microsatellite markers showed evidence for 
Hardy-Weinberg disequilibrium, allelic dropout or null 
alleles. We estimated expected heterozygosity (He) with 
GenePop on the Web (Raymond and Rousset 1995; Rous- 
set 2008) and allelic richness (Ar) using the rarefaction 
method to correct for unequal sample sizes with Fstat 3.9.3 
(Goudet 1995). We computed the population inbreeding 
coefficient (_Fis) and the average pairwise relatedness esti- 
mator (R). Fis was calculated with Genetix (Belkhir 1997) 
and reflects deviation of Hardy-Weinberg equilibrium due 
to nonrandom mating within populations. R was calculated 
with GeneAlEx 6 (Peakall and Smouse 2006) using allelic 



frequencies of the global population. Moreover, based on 
previous studies (Mendez et al. 2011), we calculated a pop- 
ulation differentiation index under a pure-drift model (_F) 
with 2MOD (Ciofi et al. 1999), which reflect the accumu- 
lated inbreeding in a population by genetic drift. Finally, 
we estimated population-specific genetic differentiation 
through the arithmetic mean of the population pairwise 
FsT (Weir and Cockerham 1984), G'^j (Hedrick 2005) and 
DsT (Jost 2008) values with Genetix and the Software for 
the Measurement of Genetic Diversity (SMOGD version 
1.2.5; Crawford 2010). 

Size and isolation of populations 

Three of the 36 surveyed localities held a small number of 
sampled individuals (n < 10). They were merged with their 
closest neighbour locality to increase sample sizes for analy- 
ses, yielding 33 local populations (Mendez et al. 2011). All 
three merged pairs of localities were in the core of the dis- 
tribution, separated from each other by 1-4 km (less than 
one-fifth of the maximum dispersal distance detected in 
the study area, Laiolo et al. 2007) and interconnected by 
suitable habitat. Moreover, pairs of merged localities did 
not show signs of genetic differentiation (low and nonsig- 
nificant pairwise _Fst) and were grouped in the same genetic 
cluster (Mendez et al. 2011). Dupont's lark population 
sizes were characterized by two different estimates: (i) the 
number of male territories (N.pop), estimated from vocal- 
izations in previous studies (Suarez 2010; Vogeli et al. 
2010), and (ii) the patch size (area) of each population. 
Area was calculated as the surface of natural steppe vegeta- 
tion occupied by the Dupont's lark, obtained from land- 
use maps and aerial orthophotographs (Vogeli et al. 2010). 
Both estimates are significantly correlated (linear regres- 
sion, = 0.616, P < 0.001, N = 33), and both are gener- 
ally good indicators of population sizes (Laiolo et al. 2008). 
They can differ, however, in small patches where the reduc- 
tion in patch size causes an increase in Dupont's lark densi- 
ties ('crowding effect'), thus impeding a stronger 
relationship between patch size and population size (Laiolo 
and Telia 2006a; Vogeh et al. 2010; Mendez et al. 2011). 
Besides, they also differ in the time and costs associated to 
data collection. N.pop estimates need exhaustive field work 
for obtaining the number of territorial males (TeUa et al. 
2005). Dupont's lark is a very elusive, secretive and difficult 
to observe terrestrial species, to the point that its occur- 
rence and distribution in Spain was not roughly assessed 
until the late 1980s (Garza and Suarez 1990). Thus, map- 
ping territorial males through their vocal activity is the only 
way to properly estimate local population sizes (Telia et al. 
2005; Laiolo and Telia 2006a). Calculating area is, however, 
relatively easy, because areas of natural steppe vegetation 
adequate for the study species are island-like remnants in a 
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landscape dominated by agriculture, so they are easily iden- 
tified with basic GIS procedures (Vogeli et al. 2010). 

We calculated population isolation by the two approxi- 
mations that were identified in a previous study (Vogeli 
et al. 2010) as the best predictors of Dupont's lark occur- 
rence at the metapopulation scale: the Euclidean distance 
fi-om each population to its nearest neighbour population 
(D.near), and the distance fi-om each population to the 
average coordinates of all Spanish Dupont's lark popula- 
tions (D. centre). While the former is straightforward to 
estimate, the latter requires a complete knowledge of the 
distribution of the species and the application of GIS tools, 
so it could be more difficult to incorporate in management 
practice. Geographical patch centres were used for calculat- 
ing both isolation indexes. 

Statistical analysis 

We used generalized linear models to assess the relation- 
ships between the population-based genetic indexes (H^, 
Ar, _Fis, R, F, F^t, Gg^, Dsx) and ecological and demo- 
graphic attributes of each local population (i.e. population 
size and isolation). Population size (characterized by N.pop 
and area) and isolation (calculated as D.centre and D.near) 
were fitted as explanatory variables. Following our hypoth- 
eses, we designed a priori five alternative models to assess 
the relationships between the explanatory variables and the 
eight response variables. The first three univariate models 
only included single explanatory variables (N.pop, area, 
D.centre or D.near) to assess their effects separately. The 
fourth model included area and isolation, while the fifth 
included area, isolation and the interaction between area 
and isolation. All models were tested twice alternatively 
using D.near and D.centre as isolation measures. We fol- 
lowed a model selection strategy using an information-the- 
oretic approach (Johnson and Omland 2004), computing 
the Akaike information criterion corrected for small sample 
sizes (AICc) and relative weight of evidence for each model 
(Wi) as the probability of model i being the best model for 
the observed data, given the set of candidate models. The 
most parsimonious model for each genetic index was 
selected based on a lower AICc and higher w,- (Johnson and 
Omland 2004). All statistical analyses were performed in R 
2.11.1 software (R Development Core Team 2011). Models 
were built with a normal distribution of errors and the 
identity link function. N.pop, area and D.near were log- 
transformed to attain normality. There was no evidence of 
overdispersion, and residuals fitted to a normal distribu- 
tion, indicating that distributions and error structures were 
appropriate (Rushton et al. 2004). We estimated parame- 
ters and plotted the finally selected, most parsimonious 
models only if the retained explanatory variables reached 
statistical significance at P < 0.05. Due to the spatial distri- 



bution of the populations, we analysed whether spatial 
autocorrelation might impact the planned analysis. We 
consequently assessed final model residuals through 
Moran's I correlograms (Dormann et al. 2007). Spatial 
autocorrelation of residuals was only significant in models 
relating H,., Fi^ and _F to N.pop and D.near, indicating they 
did not adequately capture the spatial structure of the 
genetic variance. Finally, we performed piecewise and 
exponential regression models to detect nonlinear thresh- 
olds in final models with only one independent variable 
(Toms and Lesperance 2003; Laiolo et al. 2008). Following 
AICc, piecewise models performed better in all cases, so we 
only use those for detecting thresholds. 

Results 

Our exhaustive sampling of the whole Spanish (and thus 
European) Dupont's lark distribution captured a wide 
range of both the genetic {H^ and Ar, _Fis, -R, F, _Fst> G'^j 
and Dst) and the ecological (area, N.pop, D.near and 
D.Centre) conditions (Table SI). Top ranking models for 
Ar, _R and f st and G'gj were similar to those for H^, Fis and 
Dst- For simplicity, we include the formers in Supporting 
Information and report the latters and _F in Tables 1 and 2 
and Fig. 2. Selected models included size or isolation vari- 
ables or their interaction (Table 1) and explained between 
13.6% (Dst) and 65.3% (H^) of the variation (Adj. R^) in 
the genetic parameters (Table 2). 

Population size (N.pop) was retained as a statistically 
significant predictor in final models when using D.near as 
isolation measure and explained 22% of the variance in 
heterozygosity (H^) and 24% of the variance in accumu- 
lated inbreeding (_F) (Table 2). Moreover, N.pop was also 
retained in the best competing model for Dst (14% of the 
variance explained; Table 2). Genetic diversity increased 
logarithmically, whereas genetic differentiation decreased 
logarithmically with population size (Fig. 2; panels IB, 3B). 
When analysing these relationships with piecewise regres- 
sion models, we found a significant breakpoint at 19 male 
territories in the case of He (Adj. R^ = 0.40, ^433 = 8.10, 
P < 0.001) and F (Adj. R^ = 0.25, ^4,33 = 4.52, P < 0.05) 
and at 16 for Dst (Adj. R^ = 0.21, ^433 = 3.92, P < 0.001). 
Therefore, populations with <16-19 males were prone to 
the loss of genetic diversity and to the accumulation of 
inbreeding and differentiation. 

Patch size (area) was not included as a single predictor 
but played a key role in several final models through its 
interaction with isolation (see below). 

The distance to the centroid of all populations (D.centre) 
performed better as an isolation measure than the distance 
to the nearest population (D.near) for most genetic param- 
eters (Table 2). D.centre was selected in all the genetic dif- 
ferentiation models and the accumulation of inbreeding 
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Table 1. Evaluation of alternative models for genetic parameters based on population size and isolation variables. 





Models with D. centre 








Models with D.near 










N.pop 


Area 


D. centre 


Area + 
D. centre 


Area*D. 
centre 


N.pop 


Area 


D.near 


Area + 
D.near 


Area*D. 
near 


He 

AlC, 
AAlCc 

VJj 


-91.9 
23.8 
0.0 


-85.9 
29.7 
0.0 


-99.6 
16.0 
0.0 


-103.6 
12.0 
0.0 


-115.6 
0.0 
1.0 


-91.9 
0.0 
0.6 


-85.9 
5.9 
0.0 


-87.3 
4.5 
0.1 


-89.8 
2.1 
0.2 


-89.1 
2.7 
0.1 


Fis 

AlCc 

AAIC, 

Wj 


-63.2 
15.7 
0.0 


-64.2 
14.8 
0.0 


-63.0 
16.0 
0.0 


-61.8 
17.2 
0.0 


-79.0 
0.0 
1.0 


-63.2 
6.0 
0.0 


-64.2 
5.0 
0.1 


-63.2 
6.0 
0.0 


-62.1 
7.1 
0.0 


69.2 
0.0 
0.8 


F 

AlCc 
AAlCc 

Wi 
DsT 

AlC, 
AAlCc 

Wj 


-123.0 
9.5 
0.0 

-133.2 
12.6 
0.0 


-119.5 
20.0 
0.0 

-127.4 
18.5 
0.0 


139.5 
0.0 
0.7 

145.8 
0.0 
0.7 


-137.1 
2.4 
0.2 

-143.4 
2.4 
0.2 


-134.6 
4.9 
0.1 

-141.1 
4.7 
0.1 


-123.0 
0.0 
1.0 

-133.6 
0.2 
0.4 


-119.5 
10.5 
0.0 

-127.4 
6.0 
0.0 


-121.9 
8.1 
0.0 

133.4 
0.0 
0.4 


-120.0 
10.5 
0.0 

-131.0 
2.4 
0.1 


-117.0 
13.0 
0.0 

-129.0 
4.4 
0.0 



Genetic parameters: He = expected heterozygosity, inbreeding (Fis), accumulated inbreeding (F) and differentiation (Dsj). Model selection repeated 
for two alternative measures of isolation: distance to the centroid of all populations (D. centre) and distance to the nearest population (D.near). 
N.pop = male territories, Area = patch size. Additive (+) and interactive {*) effects were considered. The corrected Akaike information criterion 
(AlCe), difference with the best model (AAICJ and the relative weight of evidence for each model (w,-) are reported (two equally good models when 
Wi < 0.5). Statistics of the most parsimonious model are highlighted in bold. 



model, explaining 25%, 40%, 41% and 43% of the variance 
of FsT, G'^j, DsT and F, respectively. On the other hand, 
D.near was selected instead of D. centre for Ar, F^i, Dsr 
and G'gj, performing similarly to N.pop (Table 1 and Sup- 
porting Information). In summary, genetic differentiation 
and the accumulation of inbreeding increased with isola- 
tion irrespectively of the isolation measure used (Fig. 2, 
panels 3 A and 4A). The piecewise regression model showed 
a significant breakpoint for Dst in 30 km to the nearest 
population (Adj. = 0.10, ^4,33 = 2.16, P < 0.001). 

Models including the interaction between patch size and 
isolation explained the highest percentages of the observed 
variance in genetic diversity (He and Ar), inbreeding (fjs) 
and relatedness (R) (65%, 55%, 42% and 24%, respectively; 
Table 2 and Supporting Information). Indeed, this interac- 
tion was retained in the most parsimonious models 
explaining inbreeding and relatedness, regardless of the iso- 
lation metric used. Figure 2 illustrates size and isolation 
interactions for (panel lA) and Fis (panel 2A,B). In all 
cases, genetic indexes varied differently with patch size 
depending on their degree of isolation, or conversely, they 
responded differently to an increase in isolation depending 
on their size. To better understand these interactions, we 
plotted the lowest and highest values of isolation found in 
our data, and the one where the effect of area did not pro- 
mote a change in genetic indexes. Above or below this level 
of isolation, the relationships between genetic indexes and 



patch area changed, either increasing or decreasing. With 
low values of isolation (23 km to the centroid or 5.6 km to 
the nearest population), genetic diversity and inbreeding or 
relatedness decreased with patch size. Distances of 103- 
145 km to the centroid of all populations or at 19 km to 
the nearest population resulted in stable genetic indexes in 
all patch sizes. At the same time, patch sizes of 300, 375, 
650 and 1800 ha made genetic indexes independent of their 
degree of isolation (intersection for R, Fis with D. centre 
and D.near and H^, respectively). Finally, larger values of 
isolation resulted in increasing genetic diversity and 
inbreeding with patch size. 

Discussion 

Our exhaustive sampling of fragmented Dupont's lark pop- 
ulations covering the whole European distribution and 
almost the 25% of the estimated male's territories, and the 
use of model selection techniques, allowed us to evaluate 
the relative importance of population (or patch) size and 
isolation in explaining the observed variance in genetic 
parameters related to diversity, inbreeding and differentia- 
tion. 

Population size had a major influence on genetic indexes 
of Dupont's larks: genetic diversity increased and differen- 
tiation decreased with local population size, being the 
expected outcome of a less intense genetic drift in larger 
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populations (Fig. 3, panel A). Nevertheless, this effect is 
likely to vary across the range of population sizes, due to 
the nonlinear relationship of these variables, and through 
time. Time is an important variable that could not be 
included directly in our analyses, but could be indirectly 
captured by one of our isolation variables (D. centre, see 
below). Dupont's lark populations holding <19 male terri- 
tories have lost genetic diversity and increase the accumu- 
lated inbreeding and above 16 have accumulated 
differentiation to detectable levels. Below these population 
size thresholds, genetic erosion was minor or undetectable 
(Fig. 3-A.3/B.3). 

The other parameter of interest, isolation, proved also of 
great importance for another genetic variable: genetic dif- 
ferentiation increased with isolation (Fig. 3, panel B). As 
predicted by theory, gene flow between habitat patches 
may buffer the effects of genetic drift. In this sense, the 
effects of isolation are expected to be stronger and more 
detectable in species with low dispersal capability (Segelb- 
acher et al. 2010) than in more dispersive species (e.g. 



Lindsay et al. 2008; Canales-Delgadillo et al. 2012), but 
also in those populations with small populations sizes suf- 
fering more intense genetic drift (Frankham et al. 2002). 
Dupont's lark presents both of these two characteristics: 
low dispersal capability and small populations sizes (Laiolo 
et al. 2007; Vogeli et al. 2010). In this study, isolation of 
local populations was measured either as the linear distance 
to the nearest population or as the distance to the centroid 
of all populations. While both measures are meant to quan- 
tify the degree of population isolation, they likely do so at 
different spatial and temporal scales. On the one hand, the 
distance to the nearest population reflects more contempo- 
rary chances in dispersal among patches. In this sense, we 
found that genetic differentiation increased faster at dis- 
tances above 30 km to the nearest population than under 
this threshold. On the other hand, the distance to the cen- 
troid of all populations performed better in reflecting long- 
term population dynamics, that is, it may better capture 
the stronger effect of genetic drift expected in edge popula- 
tions due to their longer time since isolation or to a more 
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Figure 2 Predicted relationships and partial contributions of population and ecological variables on genetic indexes: patch size {area), population size 
as number of male territories (N.pop), distance to the centroid of all Spanish populations (D.centre) and distance to the nearest population (D.near). 
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intense population size reduction and isolation. This latter 
variable is critically important, but usually missing from 
explanatory models due to the difficulty in obtaining direct 
measures (see Kyle and Strobeck 2002; Vucetich and Waite 
2003; Gapare and Aitken 2005). 

Most importantly, the interaction between patch size 
and isolation explained a much higher proportion of the 
variance in genetic diversity, inbreeding and relatedness 
than each of these variables alone, or when only additive 
effects were considered. Although the product of popula- 
tion size by isolation (captured in genetic models by the 
migration rate) is considered determinant for the spatial 
variance in gene frequency in population genetics theory, 
only a few empirical studies of wild populations have 
explicitly evaluated this interactive effect. Some studies 
considered fragmentation as the product of size and isola- 
tion, making it impossible to disentangle the relative effects 
of these interacting variables (Lange et al. 2010). Others 
included the interaction between size and isolation in the 
process of model selection, yet it was not retained during 
the model selection process (Wang et al. 2010). To our 
knowledge, this is the first study showing statistically signif- 
icant effects of this interaction on the genetic response of 
wild populations to fragmentation, and its higher explana- 
tory power with respect to the single and additive effects. 
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Interactive effects of patch size and isolation on genetic 
patterns 

As the effect of patch size on genetic patterns varied 
depending on the isolation of the patch, we discuss here 
our model predictions according to three basic scenarios of 
isolation. 

In the first scenario, where isolation is high enough to 
hamper gene flow, genetic diversity, inbreeding and related- 
ness of the Dupont's lark increased in local populations 
wdth patch size (Fig. 3, panel B). In line with our predic- 
tions, larger patches held more male territories (and thus 
more individuals) and thus retained more genetic diversity 
over time (Fig. 3-B.l). The increase in relatedness and 
inbreeding with patch size, however, is not as straightfor- 
ward to explain. Dispersion is constrained in highly isolated 
populations ( Gonzalez- Varo et al. 2010) and for our model 
species in particular, which shows extremely low propensity 
for dispersal (Laiolo et al. 2007). Hence, individuals may 
try to remain in their natal patches despite the negative 
effects for the dynamics of the population (Delgado et al. 
2011). This tendency may be strengthened by conspecific 
attraction in our study species, which may also result in 
high levels of relatedness and inbreeding, and increased 
opportunities for spatially structured physiological, genetic 
and cultural patterns within patches (Laiolo and Telia 2005, 
2008; Mendez et al. 2011; Fairhurst et al. 2013). Therefore, 
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Figure 3 General overview of the genetic consequences of patch size reduction and isolation on Dupont's lark population genetics. Based on this 
study, genetic diversity and inbreeding of Dupont's lark populations change depending of the size and isolation of steppes patches. After fragmenta- 
tion, (panel A) connected and reduced patches under 1 800 or 300 ha reduce their genetic diversity and increase their inbreeding, while (panel B) iso- 
lated (more than 13-19 km to the nearest population) and reduced patches first increase their genetic diversity and reduce their inbreeding and 
finally reduce their diversity and increase their inbreeding. We found that 16-19 males territories and 30 km to the nearest population are the frag- 
mentation thresholds needed to support the initial genetic conditions. Arrows indicate migration into or from local populations at initial, intermediate 
and final stages of patch size reduction and isolation. Shadow areas around circles indicate recent patch contractions. The number and distribution of 
birds within populations illustrate changes in population size and density throughout the stages of fragmentation. Symbols ], J. and - indicate 
increase, decrease or no change, respectively, in the corresponding genetic index: genetic diversity (He), Alelic richness (Ar), inbreeding (fis), related- 
ness (/?) or differentiation (Dst)- 



some of the increase in Fis with population size might be 
due to the sampling across genetically differentiated sub- 
populations (i.e. Wahlund effect) in larger and less dense 
patches when these are highly isolated. These populations 
might therefore favour the maintenance of overall genetic 
diversity through the formation of spatially segregated 
groups of related individuals, although inbreeding and 
average relatedness will increase. Conversely, small isolated 
populations show less genetic diversity due to the combina- 
tion of a lower effective population size and, probably, a 
longer time since isolation (Fig. 3-B.3). Surprisingly, 
inbreeding and relatedness are lower in recently contracted 
and isolated populations, which could be due to the tran- 
sient high densities occurring in these populations 
('Crowding effect'; Laiolo and Telia 2006a; Vogeli et al. 
2010; Mendez et al. 2011) favouring the contact among 
genetically differentiated groups within the patch. 



In a second scenario, we focus on populations with inter- 
mediate levels of isolation and higher levels of gene flow. 
This should theoreticaDy favour genetic diversity and 
reduce inbreeding and relatedness. In this sense, inbreeding 
and relatedness in Dupont's larks populations decreased 
with patch size (Fig. 3-A.l). Intriguingly, we detected a 
slightly increase in genetic diversity and inbreeding in well- 
connected populations with smaller patch sizes (Fig. 3-A.2/ 
A.3). This effect can be due either to a statistical artefact 
due to the relatively small sample sizes, or to the fact that 
those smaller core populations may reflect a prefragmenta- 
tion situation in which the occurrence of spatial genetic 
groups in larger patches results in heterozygote deficiency 
(due to intrapatch spatial structure), high overall diversity 
and some inbreeding (Fig. 3-B.l). 

The third scenario depicts populations where varying 
patch sizes do not cause relevant changes in genetic 
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indexes. Gene flow among these populations appears to be 
still enough to compensate for the efiiscts of drift to main- 
tain the levels of genetic diversity, and to avoid inbreeding 
and the accumulation of relatedness, irrespective of popu- 
lation size. This scenario seemed to occur in those popula- 
tions that are located <103-145 km from the geographical 
centre of the whole population or <13-19 km to the near- 
est population (for fjs and _R respectively). Remarkably, 
this value is in fair accordance with the maximum dispersal 
distance of 20 km recorded for the Dupont's lark (Vogeli 
et al. 2010). Similarly, we also identified patch sizes where 
genetic indexes are not affected by an increase in isolation. 
Approximately 300 ha are needed to maintain inbreeding 
and relatedness values irrespective of the level of isolation, 
a size which match the patch size below which Dupont's 
lark densities increase as a response to patch contraction 
(Vogeli et al. 2010). However, larger patch sizes are 
required for genetic diversity (above 1800 ha) to maintain 
levels independent of isolation. 

Conservation implications 

Some guidelines for Dupont's lark conservation can be 
extracted from the results of this study, especially regarding 
specific recommendations for minimum patch sizes and 
distances between populations needed to prevent genetic 
erosion. In terms of isolation, the distance to the nearest 
Dupont's lark population should not exceed 13 km. Popu- 
lations separated by larger distances would lose territories 
and may additionally present a crowding effect with 
numerous potentially deleterious consequences. In this 
sense, we have detected a critical distance threshold of 
30 km above which inbreeding and differentiation would 
increase dramatically. We also identified a minimum patch 
area of approximately 300 ha to avoid accumulation of 
inbreeding and relatedness, whereas larger patch sizes 
(>1800 ha, or 19 male territories) should be targeted to 
prevent the loss of genetic diversity. Although these thresh- 
olds are relative to our power to detect changes in genetic 
parameters, which is a function of sampling size, genetic 
markers and the time since fragmentation started, these 
results raise concern over the situation of the species. 
Attending to this study, the majority of Spanish popula- 
tions are currently below these threshold sizes; approxi- 
mately 90% of the occupied patches are smaller than 
1800 ha and half of them cover <300 ha (range 20- 
5000 ha). Regarding the isolation of populations, 73% of 
them are separated by more than 13 km and 54% by more 
than 19 km from the nearest population (range: 5.5- 
173 km). Consequently, many Dupont's lark populations 
may be prone to suffer an 'extinction vortex' and raise the 
possibility that genetic factors contributed to the decline in 
population's viability (Vogeli et al. 2011) and even to the 
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recent extinctions of small and isolated populations (Telia 
et al. 2005; Vogeh et al. 2010). 

Advantages and practical implications of the 
methodological approach 

The combination of ecological and genetic data in a multi- 
variate modelling framework (schematized in Fig. 4) 
allowed us to understand the effects and importance of dif- 
ferent estimators of fragmentation on population genetics 
and to detect fragmentation thresholds from which the 
genetic health of local populations would become compro- 
mised. We assessed a priori defined hypotheses on the 
effects of habitat fragmentation on genetic erosion and 
tested them through a model selection approach. Those 
hypotheses involved alternative population and ecological 
parameters and the seldom considered interactions 
between them, an aspect that can help to improve manage- 
ment strategies. The strength of our approach for under- 
standing the genetic processes derived from fragmentation 
was highlighted by the fact that finally selected models 
explained between 40% and 65% of the overall variance in 
genetic diversity, inbreeding and differentiation. Moreover, 
the use of alternative and complementary estimators of 
genetic diversity, inbreeding and differentiation (Fig. 4), 
which relate to different aspects of the genetic erosion pro- 
cesses, helped us to infer responses to fragmentation with 
different dynamics and at different spatio-temporal scales 
(Keyghobadi et al. 2005). This allows obtaining a broader 
picture of the various effects of habitat loss and fragmenta- 
tion on the genetic composition of local populations. On 
the other hand, testing several ecological variables allowed 
us to contrast and compare their association with the 
genetic indexes and make more precise management rec- 
ommendations. While previous work showed the need of 
increasing habitat size and connectivity to reverse the 
decline of this endangered species (Vogeli et al. 2010; 
Mendez et al. 2011), the fragmentation thresholds we are 
providing here offer practical guidance to wildlife manag- 
ers. 

Our conceptual and analytical framework (Fig. 4) can be 
broadened to other study systems involving any species 
threatened by habitat fragmentation. We used alternative 
variables to characterize population size and isolation, 
which are likely to differentially capture different aspects of 
the fragmentation process and that can vary between study 
systems. Number of male territories and patch size are 
alternative proxies of population size for our secretive 
study species, but they differ considerably in their estima- 
tion costs. Estimating the number of territorial males 
required costly and time-consuming censuses, while patch 
size is an easier and cheaper measure that can be obtained 
from remote-sensing or general-purpose land cover maps 
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(Fernandez 2013). However, actual population sizes, and 
even the proportion of breeding individuals, can be more 
easily obtained for more conspicuous species showing pop- 
ulation fragmentation (e.g. TeUa et al. 2013). 

Distance to the centroid of all populations and distance 
to the nearest population can be considered as complemen- 
tary measures of isolation. The former one has a 'time since 
isolation' component under our model of centripetal con- 
traction, which makes it biologically more relevant in a 
context of recent fragmentation. However, it requires the 
knowledge of the distribution of all extant local popula- 
tions, a kind of information that may not be available for 



all case studies. In this case, the use of simple and easy-to- 
measure, but stiU informative, variables (as distance to the 
nearest population) should still help to improve manage- 
ment actions. 

Finally, our conceptual and methodological approach 
can be implemented at different spatial scales. Our study 
system covered the whole European distribution of an 
endangered species, but the same approach should work 
well for studies conducted at regional or metapopulation 
scales. We thus encourage the use of our approach (Fig. 4) 
to design optimal strategies for preventing farther genetic 
erosion and for reverting it in species threatened by recent 
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and/or ongoing fragmentation. A further step to improve 
this approach would be to incorporate the explicit dynam- 
ics of habitat fragmentation and population size with time 
estimates for the onset of isolation and contraction, a type 
of information that was unfortunately not available for our 
study system. 
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centroid of all Spanish populations {D. centre) and distance to the nearest 
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